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In this paper, by using bifurcation method, we successfully find the K(2, 2) equation 
with osmosis dispersion ut + (u 2 ) x — (u 2 ) xxx = possess two new types of travelling 
wave solutions called kink-like wave solutions and antikink-like wave solutions. They 
are defined on some semifinal bounded domains and possess properties of kink waves 
and anti-kink waves. Their implicit expressions are obtained. For some concrete 
data, the graphs of the implicit functions are displayed, and the numerical simulation 
of travelling wave system is made by Maple. The results show that our theoretical 
analysis agrees with the numerical simulation. 
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1 Introduction 



In recent years, many nonlinear partial differential equations (NLPDEs) 
have been derived from physics, mechanics, engineering, biology, chemistry and 
other fields. Since exact solutions can help people know deeply the described 
process and possible applications, seeking exact solutions for NLPDEs is of 
great importance. 

In 1993, Rosenau and Hyman [1] introduced a genuinely nonlinear disper- 
sive equation, a special type of KdV equation, of the form 

u t + a(u n ) x + (u n ) xxx = 0,n> 1, (1.1) 
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where a is a constant and both the convection term (u n ) x and the disper- 
sion effect term {u n ) xxx are nonlinear. These equations arise in the process of 
understanding the role of nonlinear dispersion in the formation of structures 
like liquid drops. Rosenau and Hyman derived solutions called compactons 
for Eq.(l.l) and showed that while compactons are the essence of the focusing 
branch where a > 0, spikes, peaks, and cusps are the hallmark of the defocus- 
ing branch where a < which also supports the motion of kinks. Further, the 
negative branch, where a < 0, was found to give rise to solitary patterns having 
cusps or infinite slopes. The focusing branch and the defocusing branch repre- 
sent two different models, each leading to a different physical structure. Many 
powerful methods were applied to construct the exact solutions for Eq.(l.l), 
such as Adomain method [2] , homotopy perturbation method [3] , Exp-function 
method [4], variational iteration method [5], variational method [6,7]. In [8], 
Wazwaz studied a generalized forms of the Eq.(l.l), that is mK(n,n) equa- 
tions and defined by 

u n - x u t + a[u n ) x + b{u n ) xxx = 0, n > 1, (1.2) 

where a, b are constants. He showed how to construct compact and noncompact 
solutions for Eq.(1.2) and discussed it in higher dimensional spaces in [9]. Chen 
et al. [10] showed how to construct the general solutions and some special exact 
solutions for Eq.(1.2) in higher dimensional spatial domains. He et al. [11] 
considered the bifurcation behavior of travelling wave solutions for Eq.(1.2). 
Under different parametric conditions, smooth and non-smooth periodic wave 
solutions, solitary wave solutions and kink and anti-kink wave solutions were 
obtained. Yan [12] further extended Eq.(1.2) to be a more general form 

u m~l Ut + + b(u k ) xxx = 0,nk^ 1, (1.3) 

And using some direct ansatze, some abundant new compacton solutions, 
solitary wave solutions and periodic wave solutions of Eq.(1.3) were obtained. 
By using some transformations, Yan [13] obtained some Jacobi elliptic function 
solutions for Eq.(1.3). Biswas [14] obtained 1-soliton solution of equation with 
the generalized evolution term 

(u l ) t + a(u m )u x + b(u n ) xxx = 0, (1.4) 

where a, b are constants, while l,m and n are positive integers. Zhu et al. 
[15] applied the decomposition method and symbolic computation system to 
develop some new exact solitary wave solutions for the K(2, 2, 1) equation 

u t + (u ) x - (u ) 

xxx ~\~ u xxxxx 0, (1.5) 
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and the K(3, 3, 1) equation 

)x (li )rax "I - U X xxxx 0. 



(1.6) 



In [16], Xu and Tian introduced the osmosis K(2,2) equation 



u t + (u 2 ) x - (u 2 ) 



xxx 



0. 



(1.7) 



where the negative coefficient of dispersive term denotes the contracting dis- 
persion. They obtained the peaked solitary wave solution and the periodic 
cusp wave solution for Eq.(1.7). In this paper, we'll continue their work and 
using the bifurcation method of planar dynamical systems to derive two new 
types of bounded travelling wave solutions for Eq.(1.7). They are defined on 
some semifinal bounded domains and possess properties of kink waves and 
anti-kink waves. To our knowledge, such type of travelling wave solution has 
never been found for Eq.(1.7) in the former literature. 

The remainder of the paper is organized as follows. In Section 2, we give the 
phase portrait of the travelling wave system and use Maple to show the graphs 
of the orbits connecting with the saddle points for our purpose. In Section 3, 
we state the main results which are implicit expressions of the kink-like and 
antikink-like wave solutions. In Section 4, we give the proof of the main results. 
For some concrete data, we use Maple to display the graphs of the implicit 
functions. In Section 5, the numerical simulations of travelling wave system 
are made by Maple. A short conclusion is given in Section 6. 



2 Phase portrait of the travelling wave system 

Eq.(1.7) also takes the form 
Ut + 2uu x — Qu x u xx — 2uu xxx = 0, (2-1) 

Let u = (p(£) with £ = x — ct(c ^ 0) be the solution for Eq.(2.1), then it 
follows that 



- ciff + 2ipip' - 6ip'ip" - 2ifip f " = 0. 



(2.2) 



Integrating (2.2) once we have 



ap+(cp) 2 -2{ip') 2 -2 W " = g, 



(2.3) 
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where g is the integral constant. 



Let y = (p 



>', then we get the following planar dynamical system: 




(2.4) 



with a first integral 



H(<p, y) = ^{y 1 - -^p 2 + 3^ + = h > 



(2.5) 



where h is a constant. 

Note that (2.4) has a singular line <p — 0, to avoid the line temporarily we 
make transformation = 2ipd(. Under this transformation, Eq.(2.4) becomes 



Eq.(2.4) and Eq.(2.6) have the same first integral as (2.5). Consequently, sys- 
tem (2.4) has the same topological phase portraits as system (2.6) except for 
the straight line ip = 0. Obviously, ip = is an invariant straight-line solution 
for system (2.6). 

Now we consider the singular points of system (2.6) and their properties. 
Note that for a fixed h, (2.5) determines a set of invariant curves of (2.6). As h 
is varied, (2.5) determines different families of orbits of (2.6) having different 
dynamical behaviors. Let M(ip e ,y e ) be the coefficient matrix of the linearized 
system of (2.6) at the equilibrium point (ip e ,y e ), then 




(2.6) 




and at this equilibrium point 



we have 



J(<p e ,y e ) = det M(<p e ,y e ) 



p{<p e ,ye) = tiace(M(<p e ,y e )) 
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By the theory of planar dynamical system (see [17]), for an equilibrium point 
of a planar dynamical system, if J < 0, then this equilibrium point is a saddle 
point; it is a center point if J > and p = 0; if J = and the Poincare index 
of the equilibrium point is 0, then it is a cusp. 

Although the distribution and properties of equilibrium points of system 
(2.4) has been given in [16]. Here we also state it briefly for our purpose. 
System (2.4) has the following properties: 

(1) When g > 0, system (2.4) has two equilibrium points (<^,0) and 
(y?o,0). They are saddle points, (i) If c < 0, then there is inequality <Pq < 
| < < <p~Q ; (ii) If c > 0, then there is inequality ip^ < < | < ipQ . 

(2) When g = 0, system (2.4) has two equilibrium points (0,0) and (c, 0). 
(0, 0) is a cusp, and (c, 0) is a saddle point. 

(3) When — y < g < 0, system (2.4) has two equilibrium points (<Pq,0) 
and (<y?o,0). (i) If c < 0, then (<Po,0) is a saddle point while (y?o,0) is a 
center point. There is inequality <Pq < | < tpfi < 0; (ii) If c > , then 
(<Po,0) is a center point while (<Po,0) is a saddle point. There is inequality 

< ipo < f < </?o ■ 

2 

(4) When g = —\, system (2.4) has only one equilibrium point (|, 0). It is 
a cusp. 

(5) When g < — system (2.4) has no equilibrium point. 

Remark 2.1 Suppose that = x — ct) is a travelling wave solution for 

Eq.(1.7) for £ G (— oo, +oo) ; and lim ip(£) = A, Yon. ip{^) = B , where A and 

»— oo £-^oo 

i? are iiyo constants. If A = B, then «s called a solitary wave solution. 
If A B, then «s called a kink (or an anti-kink) solution. Usually, a 
solitary wave solution for Eq. (1.7) corresponds to a homoclinic orbit of system 
(2.6) and a periodic orbit of system (2.6) corresponds to a periodic travelling 
wave solution of Eq.(1.7). Similarly, a kink (or an anti-kink) wave solution 
of Eq.(1.7) corresponds to a heteroclinic orbit (or so-called connecting orbit) 
of system (2.6). In [16], Xu and Tian reported that when —^<g<—^-, 
system (2.6) has homoclinic orbits, and when —^<g<0, system (2.6) has 
a periodic orbit which consist of an arc and a line segment. They obtained 
peakon solutions from the limit of solitary waves and from the limit of periodic 
cusp waves. We'll obtain a new type of bounded travelling wave solutions called 
kink-like and antikink-like wave solutions for Eq.(1.7) when g > — ^§-, which 
correspond to the orbits of system (2.6) connecting with the saddle points. 

We show the phase portraits in each region and on the bifurcation curves 
in Fig.l. From Fig.l, we can see that when g > — ^j-, system (2.6) has orbits 
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connecting with the saddle points. 




(e) 



(f) 



(g) 



(h) 




(i) 



(j) 



(k) 



(0 



Fig. 1. Phase portrait of system (2.6). (a) g > 0, c < 0; (b) g > 0, c > 0; (c) 
g = 0, c < 0; (d) g = 0, c > 0; (e) < g < 0, c < 0; (f) < g < 0, 

c> 0; (g) g = , c < 0; (h) g = c > 0; (i) < 9 < -2j£, c < 0; (j) 



< 5 < 



2cf 
9 



c > 0; (k) g 



, c < 0; (1) g 



-, c> 0. 



3 Main results 



We state our main result as follows. 
Theorem 3.1 For given constant c ^ 0, let 
£ = x-ct, 



c± Vc 2 + 4# 



(3.1) 
(3.2) 



(%) When g > and c < 0, Eq.(l,7) has two kink-like wave solutions 
u = <£>i(£) and m = ^(O ^0 antikink-like wave solutions u = ^(0 and 
u = (^4(0 • 



ftfai) = /3i(a)exp(--£), Ze(-oo,&) 



(3.3) 
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PiM = A(a)exp(^), e e (-£ \°°), (3-4) 
= &(&) exp(-^), £ e oo), (3.5) 
= /3 2 (&)exp(^), £ e (-oo,e 3 ), (3.6) 



(%) VFnen g > and c > 0, Eq.(1.7) has two kink-like wave solutions 
u = <£>i(£) and n = V9 3 (£) and iwo antikink-like wave solutions u = ^(0 an< ^ 















= f3 1 (a) exp(--f), 


ee(- 




(3.7) 


AM 


= /3i(a)exp(-0, 




-e 5 ,oo)> 


(3.8) 




= /? 2 (6)exp(-i0, 




(-£o,ao), 


(3.9) 




= /? 2 (6)exp(i0, 


ee(- 




(3.10) 



When —nr < 9 < and c < ; Eq.(1.7) has a kink-like wave solution 
u = ^5(0 an ^ an antikink-like wave solution u = ^(O- 

/3 1 (^ 5 )=A(d)exp(-^), £e(-°c,£ 9 ), (3.H) 
ftfae) = A(d)exp(^), £ e (-&°o), (3.12) 



VFnen — 2|- < a < and c> 0, Eq.(1.7) has a kink-like wave solution 
u = ^7(0 and an antikink-like wave solution u = <fs(0- 

(3 2 ^ 7 )=(3 2 {k)eM-\0, ee(-oo,e o n ), (3.13) 
= exp(^), £ e (-£ n , 00), (3.14) 



where 



c — J c 2 + 3o 

^ = ^ (3-15) 



, c+Vc 2 + 3g 
V 2 = 2 


(3.16) 


h = -\{c + ^c 2 + Ag), 


(3.17) 


l 2 = l -{c 2 + Qg-c^c 2 + Ag), 


(3.18) 


h = ^(c 2 + 6g-c^c 2 + Ag), 


(3.19) 


m 1 = --(c-3y / c 2 + 4#), 


(3.20) 


m 2 = ^(c 2 + 6g + cyjc 2 + 4g), 


(3.21) 


m 3 = jj(c 2 + 6g + c^Jc 2 + 4g), 


(3.22) 


ai = c 2 + 4g -cyjc 2 + 4g, 


(3.23) 


a 2 = c 2 + Ag + c^c 2 + Ag, 


(3.24) 


6 1 = I(2c-6V / c 2 + 4 & ), 


(3.25) 


1 / 

b 2 = -(2c + Q^c 2 + 4g), 


(3.26) 


c- yfc 1 + 4g 
rvi = 

2yc 2 + 4g - c^c 2 + Ag 


(3.27) 


c + y/c 2 + Ag 
2^jc 2 + Ag + c^c 1 + Ag 




„ , \ (2vV + hip + l 2 + 2(p + h) Up - (fo ) ai 
PUW (2^^ + «iV + ^2 + &ip + / 3 ) Q1 ' 


(3.29) 


^ / x _ (2\V + m iV + m 2 + 2ip + mi)(v9 - y?o )° 2 
(2^[oo\/ lp 1 + mi 09 + mo + 6909 + m-x) 012 


(3.30) 


c 

e =2\n(p 1 (a)/p 1 (-)), 


(3.31) 


£ 3 = 21n(/3 2 (^)//3 2 (6)) 


(3.32) 


e 5 = 21n(/3 1 (a)/AK)), 


(3.33) 
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e 7 = 21n(/3 2 (|)//3 2 (6)), 



(3.34) 



e 11 = 21n(/3 2 (|)//3 2 (A;)), 



(3.35) 



(3.36) 



a, b, d, k are four constants satisfying <fi(0) = y? 2 (0) = a, ^3(0) = ^(O) = b, 
995(0) = (pe(0) = d, ^7(0) = ^8(0) = k, and there are inequalities ipQ < a < 
§ < < ip* 2 < b < for c < 0, (Pq < a < </?* < < § < b < ip$ for c > 0, 
(Pq < d < § < y?o < for c < and < ip§ < ~ < k < ip£ for c> 0. 

We will give the proof of this theorem in Section 4. Now we take a set of data 
and employ Maple to display the graphs of u — <pi(£)(i — 1, 2, 3, 4, 5, 6, 7, 8). 

Example 3.1 Taking g = 5, c = —1 (corresponding to (1) in Theorem (3.1)), 
it follows that (fo = -2.79129 , tp$ = 1-79129, h = -4.24924, l 2 = 5.93043, 
l 3 = 23.7217, ai = 25.5826, 61 = -9.83182 , a± = -0.551865. Further, choos- 
ing a = —0.75 G (ipo, f ), we obtain Q = 0.0482492. We present the graphs of 
the solutions <fi(£) and <^ 2 (0 i n Fig. 2 (a) and (b), respectively. Meanwhile, we 
get mi = 4.51591, m 2 = 4.4029, m 3 = 17.6116, a 2 = 16.4174, b 2 = 8.49848, 
a 2 = 0.442092, <p\ = 1.5. Further, choosing b = 1.6 G (ip2,<Po), we get 
= 0.343656. The graphs of the solutions ^(0 and (f^) are presented in 
Fig. 2(c) and (d), respentively. The graphs in Fig. 2 show that <p\(£) and y?3(£) 
are two kink-like wave solutions and <^ 2 (£) and are two antikink-like 

wave solutions. 

Example 3.2 Taking g = 5, c = 1 (corresponding to (2) in Theorem (3.1)), it 
follows that v?o = -1-79129 , tp$ = 2.79129, h = -4.51591, l 2 = 4.4029, l 3 = 
17.6116, ai = 16.4174, h = -8.49848 , on = -0.442092, ip* = -1.5. Further, 
choosing a = — 1.6 G (ipo, <£*), we obtain £q = 0.343656. We present the graphs 
of the solutions <pi(£) and y? 2 (£) in Fig. 3 (a) and (b), respectively. Mean- 
while, we get m x = 4.249241, m 2 = 5.93043, m 3 = 23.7214, a 2 = 25.58264, 
b 2 = 9.93182, a 2 = 0.551865. Further, choosing b = 2 G (|,y?o), we get 
£0 = 0.773847. The graphs of the solutions <ps(£) and <^4(0 o,re presented in 
Fig. 3(c) and (d), respentively. The graphs in Fig. 3 show that <pi(£) and <ps(£) 
are two kink-like wave solutions and </? 2 (£) and are two antikink-like 

wave solutions. 

Example 3.3 Taking g = —0.5, c = —2 (corresponding to (3) in Theorem 
(3.1)), it follows that ifo = -1.70711 , (ft = -0.292893, h = -0.747547, 
l 2 = 0.638071, l 3 = 2.5528, a 1 = 4.82843, h = -4.16176 , a x = -0.776887. 
Further, choosing d = —1.2 G (</?o > §); we obtain £q = 0.448123. We present 
the graphs of the solutions <^s(0 and <pe(£) * n Fig -4 ( a ) an d (b), respectively. 
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0.5 q Oil 



-7 "6 -5 -4 -3 -2 -1 ff 12 3 4 5 6 7 

5 -0.5^, -0.5! q 




T7ie graphs in Fig. 4 show that <£>s(£) a kink-like wave solution and <fe(0 is 
an antikink-like wave solution. 



Example 3.4 Taking g = —0.5, c = 2 (corresponding to (4) in Theorem 
(3.1)), it follows that <po = 0.292893 , <p£ = 1.70711, m x = 0.747547, m 2 = 
0638071, m 3 = 2.55228, a 2 = 4.82843, b 2 = 4.1676, a 2 = 0.776887. Further, 
choosing k — 1.2 £ (fj^o); u ' e 9 e ^ ^o 1 = 0.448123. T/je graphs of the solutions 
<P7(£) and (fis{0 are presented in Fig. 5 (a) and (b), respentively. The graphs in 
Fig. 5 show that (fj^) is a kink-like wave solutions and (ps(£) is an antikink-like 
wave solutions. 



4 Proof of main results 



Suppose g > and c < 0, then system (2.4) has two saddle points (<p , 0) 
and (tpo, 0). There are four orbits connecting with (<Pq, 0). We use I 2 - to denote 

the two orbits lying on the right side of ((po , 0) (see Fig.6(a)). Meanwhile, there 
are four orbits connecting with (<Pq, 0). We employ l 1 + and l 2 + to denote the 

two orbits lying on the left side of (<Po,0) (see Fig.6(a)). 
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(a) 



1 2 3 4 5 



-1 

-1.60 
h -1.65 
-1.70- 
-1.75- 



(b) 



1 2 3 4 5 




(c) (d) 
Fig. 3. The graphs of tpi(£)(i = 1, 2, 3, 4) when g = 5, c = 1, a = —1.6, b = 2. 



-1 -6 -5 -4 -3 -2 -1 





(b) 



Fig. 4. The graphs of <Pi(£)(i = 5,6) when g = —0.5, c = —2, <i = —1.2. 

3 7l ;2 ;i ov1 ^ 72 

expressions, respectively 



On the oj — y plane, the orbits Z 1 -, Z 2 _, /* + and I 2 + have the following 

V? <Po 



z 1 - : y 



^0 ' ^ 



((fo - v)Vv' 2 + h<P + k 
(if - y?Q )^/(f 2 + hip + / 2 

209 



(4.1) 
(4.2) 
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(c) (d) 

Fig. 6. The sketches of orbits connecting with saddle points. (a)g > 0, c < 0; 
2 4<9<0, c<0; (d)-^ 



(b)g > 0, c> 0; (c) -\ < g < 0, c < 0; (d)-%- < 5 < 0, c> 0. 



,i (<Po-vW<P 2 + miy? + m7 , , 

' ^ 2y? 1 ' J 



y 2<^ v y 



where y2 and are in (3.2) , l-y and ^ are i n (3-17) and (3.18), mi and wi2 
are in (3.20) and (3.21), respectively. 
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Assume that <p = <£i(£)> V 9 = ^(O) V 9 = ^3(0 an d V 9 = ^(O 011 

^0 

/ 2 _, / x + and l 2 + , respectively and y?i(0) = ^2(0) = a, ^(O) = y?4(0) = b, 

^0 ^0 fo 

(^2 = C+ ^2 +3g '; where a and 6 are two constants satisfying ip^ < a < | and 
<P2 < b < <£q. Substituting (4.1)-(4.4) into the first equation of (2.4) and 
integrating along the corresponding orbits, respectively, we have 

/? 7L; ds=-fds (along £), (4.5) 

a 

[- s * =ds = - ds (along I 2 .), (4.6) 

6 -s 1 ° 

/ 7 + w 2I ^=ds = - [ds (along f + (4.7) 



/ 7 +\ 1 L ■ = d3 = 9 / ds ( alon s ■ ( 4 - 8 ^ 

J (S — (^oJV s + ?TliS + 7712 2 j ^0 



With the aim of Maple, we obtain the implicit expressions of <fi(£) as in (3.3)- 
(3.6). 

Meanwhile, suppose that <£>i(£) — » f as £ — > £q, ^(O — * f as £ — * ~~ £o> 
^(O — ^ 9^2 as C — ^ ~~ £0 an d ^4(0 ~~ > ¥2 as £ — * £0' t nen ^ follow from (4.5) 
-(4.8) that 



2 

£0 = Co = / 7 rr 7' . =, ds (along |J ), (4.9) 

J (s - Po)v s + hs + h 



V>2 

Co 3 = ^ = /" 7 , = = ds (along (4.10) 
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With the aim of Maple, we get the expressions of £g and £g as i n (3.31) and 
(3.32). The proof of (1) in Theorem (3.1) is finished. 

Similarly, we can prove (2)-(4) in Theorem (3.1). Here we omit the details. 
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5 Numerical simulations 



In this section, we will simulate the planar graphs of the kink-like and 
antikink-like wave solutions. 

From Section 2, we see that in the parameter expressions ip = </?(£) and 
V = y(C) °f t ne orbits of system (2.6), the graph of </?(£) and the integral 
curve of Eq.(2.3) are the same. In other words, the integral curves of Eq.(2.3) 
are the planar graphs of the traveling waves of Eq.(1.7). Therefore, we can 
see the planar graphs of the kink-like and the antikink-like waves through the 
simulation of the integral curves of Eq.(2.3). 

Example 5.1 Take the same data as Example (3.1), that is g = 5, c = — 1, 
a = —0.75, b = 1.6. Let if = a = —0.75 in (4-1) an d (4-%)> then we can get 
y w 4.23397 or y w -4.23397. And let ip = b = 1.6 in (4-3) and (44), then 
we obtain y pa 0.23198 or y pa —0.23198. Thus we take the initial conditions 
of Eq.(2.3) as follows: (a) Corresponding to we take <p(0) = —0.75 and 

<p'(0) = 4.23397; (b) Corresponding to / 2 _ we take tp(Q) = -0.75 and <p'(0) = 
-4.23397; (c) Corresponding to we take ip(0) = 1.6 and <p'(0) = 0.23198; 
(d) Corresponding to Z 2 + we take <p(0) = 1.6 and <p'(0) = —0.23198. Under 
each set of initial conditions we use Maple to simulate the integrals curve of 
Eq.(2.3) m Fig.7. 

Example 5.2 Take the same data as Example (3.2), that is g = 5, c = 1, 
a = —1.6, 6 = 2. Let ip = a = —1.6 in (4-1) and (4-2) , then we can get 
y pa 0.230189 or y pa -0.230189. And let ip = b = 2 in (4.3) and (4.4), then 
we obtain y pa 0.849228 or y pa —0.849228. Thus we take the initial conditions 
of Eq.(2.3) as follows: (a) Corresponding to r_ we take y?(0) = —1.6 and 

ip'(0) = 0.230189; (b) Corresponding to Z 2 _ we take <p(0) = -1.6 and ip'(0) = 

^0 

-0.230189; (c) Corresponding to I 1 + we take <p(0) = 2 and ip'(0) = 0.849228; 
(d) Corresponding to / 2 + we take <p(0) = 2 and <p'(0) = —0.849228. Under 
each set of initial conditions we use Maple to simulate the integrals curve of 
Eq.(2.3) in Fig. 8. 

Example 5.3 Take the same data as Example (3.3), that is g = —0.5, c = 
—2, d = —1.2. Let <p = d = —1.2 in (4-1) and (4-2), then we can get y pa 
0.364453 or y pa —0.364453. Thus we take the initial conditions of Eq.(2.3) as 
follows: (a) Corresponding to I 1 _ we take <p(0) = —1.2 and (f'{0) = 0.364453; 

(b) Corresponding to Z 2 _ we take <p(0) = —1.2 and <p'(0) = —0.364453. Under 
each set of initial conditions we use Maple to simulate the integrals curve of 
Eq.(2.3) in Fig. 9. 
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phi[3](xi) 



phi[4](xi) 



-1 1 2 3 4 5 



-6 -5 -4 -3 -2-10 1 



(c) 



(d) 



Fig. 7. The numerical simulations of integral curves of Eq.(2.3) when g = 5 and 
c = -1. (a) <p(0) = -0.75, <p'(0) = 4.23397; (b) 99(0) = -0.75, y>'(0) = -4.23397; 
(c) 99(0) = 1.6, </>'(0) = 0.230189; (d) 99(0) = 1.6, tp'(0) = -0.230189. 

Example 5.4 Take the same data as Example (3.4), that is g = —0.5, c = 
2, k = 1.2. Lei ip = k = —1.6 m and (C/.^J , i/ien we can get y m 

0.364453 or y ~ —0.364453. Tans we £a£;e tae initial conditions of Eq.(2.3) as 

1.2 and p'(O) = 0.364453; 

we ta£;e y?(0) = 1.2 and yj'(0) = -0.364453. Under 
each set of initial conditions we use Maple to simulate the integrals curve of 
Eq.(2.3) in Fig. 10. 



follows: (a) Corresponding to we take <f(0) 
(b) Corresponding to 



Comparing Fig. 2 with Fig. 7, Fig. 3 with Fig. 8, Fig.4 with Fig.9, and Fig. 5 
with Fig. 10, we can see that the graphs of <fi(£)(i = 1, 2, 3, 4, 5, 6, 7, 8) are the 
same as the corresponding integral curves of Eq.(2.3). This implies that our 
theoretical results agree with the numerical simulations. 



6 Conclusion 



In this paper, we find a new type of bounded travelling wave solutions 
for the K(2, 2) equation with osmosis dispersion. Their implicit expressions 
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xi 

1 2 3 4 5 6 




(c) (d) 

Fig. 8. The numerical simulations of integral curves of Eq.(2.3) when g = 5 and 
c = 1. (a) p(0) = -1.6, <p'(0) = 0.230189; (b) </?(0) = -1.6, ip'(0) = -0.230189; (c) 
(p(0) = 2, ip'(0) = 0.849228; (d) (^(0) = 2, ^'(0) = -0.849228. 




are obtained in (3.3)-(3.14). From the graphs (see Fig.2-Fig.5) of the implicit 
functions and the numerical simulations (see Fig.7-Fig.10) we see that these 
new bounded solutions are defined on some semifinal bounded domains and 
possess properties of kink and anti-kink wave solutions. 



1(3 




(a) 




(b) 



Fig. 10. The numerical simulations of integral curves of Eq.(2.3) when g = —0.5 
and c = 2. (a) 99(0) = 1.2, p'(0) = 0.364453; (b) ip{0) = 1.2, ip'(0) = -0.364453. 
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